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Cn| \ Abstract. Adopting the definition of excedances of type B due to Brenti, 

Q-ij we give a type B analogue of the g-derangement polynomials. The connec- 

tyy • tion between g-derangement polynomials and Eulerian polynomials naturally 

extends to the type B case. Based on this relation, we derive some basic prop- 
erties of the g-derangement polynomials of type B, including the generating 
q | function formula, the Sturm sequence property, and the asymptotic normal 

distribution. We also show that the g-derangement polynomials are almost 
symmetric in the sense that the coefficients possess the spiral property. 
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^ : 1 Introduction 

q 

In this paper, we define a type B analogue of the g-derangement polyno- 
00 ! mials introduced by Brenti [1] by g-counting derangements with respect to 

the number of excedances of type B, also introduced by Brenti [3]. We give 
some basic properties of these polynomials. It turns out that the connection 
between the g-derangement polynomials and the Eulerian polynomials nat- 
urally extends to the type B case, where the type B analogue of Eulerian 
polynomial has been given by Brenti [3], and has been further studied by 
Chow and Gessel in [6]. 



Let us now recall some definitions. Let & n be the set of permutations of 
[n] = {1,2,..., n}. For each a G &„., the descent set and the excedance set 
of a = (J\<Ji ■ • • o n are defined respectively as follows, 

Des(cr) = {% G [n - 1] : cr* > cr m }, 

Exc(cr) = {i G [fi - 1] : a, > i}. 

The descent number and excedance number are defined respectively by 

des(er) = |Des(a)|, exc(cr) = |Exc(c)|. 
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The Eulerian polynomials [TTI [13] are defined by 

A n (q) = q dcs(a)+1 = J2 Q CXC(a)+ \ n>0. 



The Eulerian polynomials have the following generating function 

t n (1 - q)e' 

!-r = — — 

n! e qt — qe 



^ . . .t n {l-q)e qt 



n>0 



A permutation a = o"iO"2 • ■ - <7„ is a derangement if <jj 7^ z for any 16 [n]. 
The set of derangements on [n] is denoted by D n . Brenti [1] defined the 
g-derangement polynomials of type A by 



d n (q) = } , q 



^ ,exc(cx) 



ctGD„ 



and proved that <i„(g) is symmetric and unimodal for n > 1. The following 
formulas (11. 2p and (II. 3p are derived by Brenti [I]. 

Theorem 1.1 For n > 0, 

4(g) = £(-ir*(™W). (1.2) 

fc=0 ^ ' 



Theorem 1.2 PFe /jai>e 

E d ^-\ = 1 _ En>2 ( g + g 2 + . . . + q n-l )t n/ n \ ■ (°) 



n>0 



A combinatorial proof of the above identity is given by Kim and Zeng 
[S] based on a decomposition of derangements. Brenti further proposed the 
conjecture that d n (q) have only real roots for n > 1, which has been proved 
independently by Zhang [TJ], and Canfield as mentioned in [2]. 

Theorem 1.3 The polynomials d n (q) form a Sturm sequence. Precisely, 
d n (q) has n distinct non-positive real roots, separated by the roots of d n -i(q). 



The following recurrence relation is given by Zhang |14j . which has been 
used to prove Theorem 11.31 

Theorem 1.4 For n > 2, we have 

dn{q) = (n - l)qd n -i(q) + q(l - q)d' n _^{q) + (n - i)qd n - 2 (q). 
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This paper is motivated by finding the right definition of a type B ana- 
logue of the g-derangement polynomials of type A so that we can get anal- 
ogous properties to the above theorems for the type A g-derangement poly- 
nomials. We discover that the notion of excedances of type B introduced 
by Brenti serves as the right choice for type B derangement polynomials, 
although there are several possibilities to define type B excedances, see 
[31 [12] • Nevertheless, it should be noted that the type B derangement 
polynomials are not symmetric compared with the case of type A. On the 
other hand, we will be able to show that they are almost symmetric in the 
sense that their coefficients have the spiral property. 

This paper is organized as follows. In Section 2, we recall Brenti's def- 
inition of type B excedances, and present the definition of g-derangement 
polynomials of type B, denoted by d^(q). In Section 3, we establish the 
connection between the derangement polynomials of type B and the Eule- 
rian polynomials of type B. This leads to a generating function formula for 
type B derangement polynomials. We then extend the [/-algorithm and V- 
algorithm given by Kim and Zeng [8] to derangements of type B. This gives 
a combinatorial interpretation of the generating function formula. In Section 
4, we prove that the polynomials d^(q) form a Sturm sequence. We also 
show that the coefficients of d„(q) possess the spiral property. In Section 5, 
by using Lyapunov's theorem we deduce that the limiting distribution of the 
coefficients of d„(q) is normal. 

2 The Excedances of Type B 

In this section, we recall Brenti's definition of type B excedances, and give 
the definition of the g-derangement polynomials of type B. We adopt the 
notation and terminology on permutations of type B, or signed permutations, 
as given in [5]. Let B n be the hyperoctahedral group on [n]. We can view 
the elements of B n as signed permutations of [n], written as a = a±a 2 ■ ■ ■ o~ n , 
in which some elements are associated with a minus sign. We also express a 
negative element — i in the form i. 

The type B descent set and the type B ascent set of a signed permutation 
a are defined by 

Des B (cr) = {i G [0, n - 1] : o* > a i+ i}, 
AscbO) = {i e [0,n - 1] : <7i < a i+1 }, 

where we set o"o = 0. The type B descent and ascent numbers are given by 

dese(cr) = |Des£(cx)|, ascs(cr) = |Asc#(er)|. 
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A derangement of type B on [n] is a signed permutation a = o\Oi • • • a n 
such that o~i ^ i, for all i G [n\. The fixed point is a position i such that 
o~i = i. The set of derangements in B n is denoted by D%. 

Let us recall the definitions of excedances and weak excedances of type 
B introduced by Brenti [3]. For further information on statistics on signed 
permutations, see [3JE1E]. 

Definition 2.1 Given a G B n and i G [n], we say that i is a type B ex- 
cedance of a if either o~i = —i or o~\ ai \ > o~i. We denote by excs(cr) the 
number of type B excedances of a. Similarly, we say that i is a type B weak 
excedance of a if either o~i — i or a\ Ci \ > o~i, and we denote by wexcs(cx) the 
number of type B weak excedances of a. 

Based on the above definition of type B excedances, we define the type 
B analogue of the g-derangement polynomials. 

Definition 2.2 The type B derangement polynomials d B (q) are defined by 

n 

d B n {q) = <? CXCB(CT) = E<^*> n ^ X > 

where d n> & is the number of derangements in D B with exactly k excedances 
of type B. For n = 0, we define d^(q) = 1. 

are the polynomials d^(q) for n < 10: 

q- 

Aq + g 2 , 
8q + 20q 2 + q 3 , 
16q + lUq 2 + 72q 3 + q A , 
32q + 752g 2 + 1312g 3 + 232g 4 + q 5 , 
64g + 3456g 2 + 14576g 3 + 9136g 4 + 716g 5 + q 6 , 
128g + 14912g 2 + 127584g 3 + 190864g 4 + 55624g 5 + 2172g 6 + q\ 
256g + 62208g 2 + 977920g 3 + 2879232g 4 + 2020192g 5 
+ 314208g 6 + 6544g 7 + g 8 , 

512g + 254720g 2 + 6914816g 3 + 35832320g 4 + 49168832g 5 
+ 18801824g 6 + 1697408g 7 + 19664g 8 + q 9 , 



Below 

df(q) = 

<${q) = 

di(q) = 

d B Aq) = 

d B {q) = 

d B (q) = 

d B (q) = 

di(q) = 

d B (q) = 
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df (q) = 1024g + 1032192g 2 + 46429440g 3 + 394153728g 4 + 937670016g 5 
+ 704504832g 6 + 161032224g 7 + 8919456g 8 + 59028g 9 + g 10 . 



3 The Generating Function 

The first result of this section is a formula expressing d^(q) in terms of 
B n (q), where B n (q) are the Eulerian polynomials of type B. This formula 
is analogous to that of Brenti [1] for the type A case , and it enables us to 
derive a formula for the generating function of d^{<i)- 

The Eulerian polynomials of type B are defined by Brenti [3] based on 
the number of descents of type B: 



Brenti [3] obtained the following formula for the generating function of the 
Eulerian polynomials of type B, see, also, Chow and Gessel [B], 



The following theorem is obtained by Brenti [3], which will be used to 
establish the formula for d^(q). 

Theorem 3.1 There is a bijection ip: B n — > B n such that 



for any a G B n . 

The following formula indicates that the notion of excedances of type B 
introduced by Brenti is a right choice for type B derangement polynomials. 

Theorem 3.2 We have 




(3.1) 




(3.2) 



ascg (<^(cr)) = wexcs(cr) 




k=0 



(3.3) 
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Proof. From Theorem I3.1[ we see that the number of excedances of type B 
and the number of descents of type B are equidistributed on B n . So we 
deduce that 



It suffices to construct a correspondence between B n and {(S, T, a)\S U T = 
[n] , S R T = 0and<r G B| S |} such that excs(7r) = excs(cr) for any tt G £> n . 

Given a signed permutation, we can decompose it into two parts separat- 
ing the fixed points from the non-fixed points. Precisely, each ir G B n can be 
represented by (S, T, a), where S is the set of non-fixed points of ir, T is the 
set of fixed points of ir, and a is a reduced signed derangement of ir. Keep 
in mind that i is a fixed point of ir if 7Tj = i. Let \S\ — k, then a is obtained 
from 7r as a signed derangement on [k] by deleting the fixed elements and 
reducing the resulting signed permutation to [k] . Formally speaking, let r be 
the signed permutation obtained from it by deleting the fixed points, then a 
is derived from r by replacing the minimum element by 1 or 1 depending on 
its sign, and replacing the second minimum element by 2 or 2 depending on 
its sign, and so on. Note that the elements in r are ordered regardless of their 
signs. For example, let it = 4 6 3 75 1 2. Then S = {1, 2, 4, 6, 7}, T = {3, 5}, 
r = 467l2and(T = 345l2. 

On the other hand, let (S, T, a) be a representation of a signed permuta- 
tion on [n] such that S U T = [n], S fl T = 0, |5| = and cr is a signed de- 
rangement on [k]. Then we can recover a unique signed permutation ir on [n]. 
Let S = {si, s 2 , . . . , s fc } with s 1 < s 2 < ■ ■ ■ < s k and T = {t 1} t 2 , . . . , t„_ fc } 
with ti < t 2 < ■ ■ ■ < t n -fc- First r can be obtained from a by replacing % or 
i in a by Sj or Sj, then 7r can be constructed from r by inserting the fixed 
points. 

We proceed to show that excs(7r) = excs(cr). Assume that % is an ex- 
cedance of ir, that is, tti = —i or Tr^i > 7r». Clearly, the fixed points do 
not contribute to the number of excedances. Then there exists j such that 
7Tj = Tj. If 7Tj = — i = Tj, then there are i — j fixed points in ir before i. 
This implies that Tj = —i is the j-th minimum element in r regardless of 
the signs. By the transformation from r to o as given before, we see that 
<Tj = — j, that is, j is an excedance of a. 

We now come to the case > 7Tj. Clearly, i is not a fixed point; oth- 
erwise, we would deduce that = 7r, ; > ?Tj. We have two cases depending 




(3.4) 



We will establish the following relation 




(3.5) 
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on the sign of 7Tj. First, we assume that 7Tj is positive. Then 7Tj is not a fixed 
point of 7T. For notational convenience, let 7Tj = j, so we have iij > j. Since 
neither i nor j is a fixed point of 7r, both should appear in r. Thus there 
exist ji and ji such that 7Tj = Tj l and 7Tj = Tj 2 . With the above notation, we 
see that Tj 2 > Tj 1 . This implies that <jj 2 > <jj l , since the transformation from 
t to o~ is order preserving. In view of 7Tj = j and ttj = Tj 2 , we find that 7Tj = j 
is the j2-th minimum element in r regardless of the signs, namely, aj 1 = j 2 . 
Hence we deduce that a aj = aj 2 > Oj x , namely, ji is an excedance of o with 
Gj 1 being positive. Conversely, given any excedance ji of a with <Tj l being 
positive, we can reverse the above procedure to generate an excedance i of 7r 
with 7Tj being positive. 

It remains to consider the case when 7Tj = j is negative. It is clear that i 
is not a fixed point of 7r. With this notation, we have ttj > j. Since i is a not 
a fixed point of 7r and 7Tj is negative, both i and j will appear in r. Hence 
there exist ji and j'2 such that 7Tj = and 7Tj = Tj 2 . Moreover, we see that 
Tj 2 > Tj 1 and <7j 2 > Oj x . Using the same procedure as given before, we find 
that <jj x = j2- It follows that cr\ (Tj \ = Gj 2 > er^, i.e., ji is an excedance of 
a with Oj x being negative. Conversely, given an excedance ji of a with 
being negative, we can reverse the above procedure to generate an excedance 
i of 7r with 7Tj being negative. 

Combining the above cases, we arrive at the conclusion that excs(7r) = 
exc#(cr), which implies (13.51) . Hence we get the following relation 

B n(^EQft). (3-6) 

Using the binomial inversion, we arrive at (13. 3p . This completes the proof. I 

The generating function of d^(q) is then obtained from the generating 
function of B n (q). 

Theorem 3.3 We have 

Td B (a) tn - ^-^ - e - (37) 



n>0 



Proof. Using (13.21) and (I3.6p . we obtain 

, R/ J n . J n (l-q)e t( - 1 -^ , . 

n>0 ' n>0 * 

The last equality of (13.71) is straightforward. This completes the proof. I 

Next, we give a combinatorial interpretation of the identity (13. 7p based 
on a generalization of the decomposition of derangements given by Kim and 
Zeng [S] for their combinatorial proof of (11.31) . 
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A Combinatorial Proof of Theorem \3.3\ . First, we give an outline the proof 
of Kim and Zeng for ordinary derangements. We adopt the convention that 
a cycle a = sis 2 ■ ■ ■ Sk of length k is written in such a way that s± is the 
minimum element and a Si = s i+ i with s fc+1 = si- A cycle a (of length at 
least two) is called unimodal (resp. prime) if there exists i (2 < % < k) 
such that Si < ■ ■ ■ < Sj_i < Sj > s i+ i > ■ ■ ■ > Sk (resp. in addition, 
< Sk). Let (li, . . . , l m ) be a composition of n, a sequence of prime cycles 
r = (n, r 2 , . . . , r TO ) is called a P-decomposition of type (/ 1; . . . , l m ) if 7$ is 
of length li and the underlying sets of Ti,r 2 , . . . ,r m form a partition of [n]. 
Define the excedance of r as the sum of the excedances of its prime cycles, 
that is, 

exc(r) = exc(ri) + • • • + exc(r m ), 

and the weight of r is defined by q exc ( T \ Note that one needs to express 
a cycle in the two row permutation form for the purpose of computing the 
excedances. The details are omitted here. In [8], Kim and Zeng obtained 
a bijection which maps the number of excedances of a derangement to the 
number of excedances of a P-decomposition with type (h, . . . , l m ), U > 2. 
Then the generating function of d n (q) follows from the generating function 
of P-decomposition with type (/1, . . . , l m ), as given by 

//, J |_ / \ JUL fh+-+lm 

Summing over h, ■ ■ ■ ,l m > 2 and m > 0, we are led to the right hand side of 
the generating function formula (11. 3j) . 

We now proceed to extend the above construction to type B derange- 
ments. We need the cycle decomposition of a signed permutation, which can 
be viewed as the cycle decomposition of an ordinary permutation with signs 
attached to some elements. There is one point that needs to be taken into 
account, that is, a signed permutation is a signed derangement if and only 
if the cycle decomposition does not have any one-cycle with a positive sign. 
More precisely, for any derangement tt of type B, we can decompose it into 
cycles 

71 = (Ci, C 2 , . . . , Cfc), 

where Ci,C 2 , ■ ■ ■ ,Ck are written in decreasing order of their minimum ele- 
ments subject to the following order 

n<---<2<I<l<2<---<n. (3.9) 

We first apply the [/-algorithm [S] to decompose each cycle a = Sis 2 ■■ ■ Sk 
of n into a sequence of unimodal cycles, here we impose the order (13.91) in 
defining unimodal cycles. It should be noted that a cycle with only one 
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negative element is also considered as a unimodal and prime cycle. Then we 
define U(n) = (17(d), U(C 2 ), . . . , U(C k )). 

The [/-algorithm 

1. If a is unimodal, then set U(a) = (a). 

2. Otherwise, let i be the largest integer such that Sj_i > Sj < Sj + i and 
j be the unique integer greater than % such that Sj > Sj > Then 
set £7(er) = (£/(cri), er 2 ), where a\ = S\ • • • Sj_iS J+ i • • • s^, and cr 2 = 
SjS i+ i • • • Sj is unimodal. 

We claim that the number of excedances of tt is equal to total number 
of excedances of unimodal cycles in U(tt). It suffices to verify that this 
statement is valid for each cycle a = Sis 2 ■ ■ ■ of 71". Clearly, it is true if a 
is unimodal. Otherwise, it suffices to show excs(cx) = excs(o"i) + excs(cr 2 ). 
Assume \s t \ is an excedance of a, i.e., cr| CT| | > C| St |. Then it is necessary to 
find an excedance in o\ or er 2 . By the cycle notation of a, we have a\ St \ = s t +i 
and <T\ a , ,| = s 4+ 2, then |s t | is an excedance of a implies s t +2 > St+i- For 
t — i — 2, we have s t+2 = s« and s t+ i = Sj_i. On the other hand, we have 
Sj < by the choice of 2, so it cannot be an excedance of a. Using the 
same argument, we see that when t — j — 1, \st\ cannot be an excedance of 
a. Therefore, ifl <t<i — 2orj + l <t<k, then \st\ is an excedance of o\. 
Similarly, if % < t < j — 1, then \s t \ is an excedance of <r 2 . If t = i — i.e., 
is an excedance of a, then Si + \ > Si, which implies that \sj\ is an excedance 
of cr 2 . If t — j, i.e., \sj\ is an excedance of a, then s-,- +2 > Sj+i, which implies 
that is an excedance of a\. Conversely, given an excedance of o\ or cr 2 , 

we can determine an excedance of a by reversing the above procedure. 

For example, let 7r = 354296871. Then we have excs(vr) = 5, and 
d = 7 8, d = 6 and C 3 = 5 9 1 3 4 2. Moreover, we find 

U{C X ) = (78), U(C 2 ) = (6), f/(d) = (59, 1342), 

and 

U(tt) = (78, 6, 59, 1342). 

Note that excB{U(7r)) = 5, in accordance with exc£(7r) = 5. 

Next, we recall the U-algorithm given in [8], which transforms a sequence 
of unimodal cycles into a sequence of prime cycles. For signed derangements, 
we will use this algorithm by imposing the order relation (13.91) . 

The V-algorithm 

1. If cr is prime, then set V(a) = (a). 
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2. Otherwise, let j be the smallest integer such that Sj > Sj > Sj+% > Sj_i 
for some integer i greater than 1. Then set V(a) = (V(ai), 02), where 
o\ = s\ ■ • • Si-iSj + \ ■ ■ ■ Sk, and a 2 = SjSj+i • • • Sj is prime. 

We claim that the total number of excedances of unimodal cycles in U(ir) 
is equal to total number of excedances of prime cycles in V o U (n) . It suffices 
to prove that the claim is valid for any unimodal cycle a = S1S2 ■ ■ ■ s&. Clearly, 
it is true if a is prime. Otherwise, without loss of generality we may show 
that excs(cx) = exc^Ci) + exc^o^). Assume \s t \ is an excedance of a, we 
have s t+ 2 > s t . As will be shown, we can find an excedance of a% or a 2 . For 
t = j — 1, | St | cannot be an excedance of a, since Sj + i < Sj by the choice of 
j. So, if 1 < t < z — 2oij + l<t<k, then \st\ is an excedance of o\. If 
i < t < j — 1, we deduce that \s t \ is an excedance of Since Si > Sj_i by 
the choice of i, when t = i — 2, \st\ is an excedance of a. On the other hand, 
we have Sj+x > by the choices of % and j, which implies that \st\ is an 
excedance of <j\ . For t = i — 1 or t = j , we can use the same argument as in 
the [/-algorithm. 

Applying V to each cycle of U(tt) in the above example, we obtain that 
V oU(n) = (78, 6, 5 9, 12, 3 4), 

and excsiV o U(tt)) = 5. 

Using the composition V o U, we can transform a derangement in B n 
to a P-decomposition of [n]. Moreover, it has been shown that this map 
is a bijection in [H]. In the type B case, we define the weight of each 
prime cycle r by g oxc s(^) ; w here excetV) is the number of the excedances 
of the type B derangement with cycle decomposition r. Note that in the 
cycle decomposition of a type B derangement, we allow cycles of length one 
with negative elements. Thus the corresponding P-decompositions have type 
(l fc , h, . . . , l m ), k > 0, U > 2. For a cycle containing only one negative ele- 
ment, the weight is q. For a cycle of length I > 2, we have 2 l choices for the I 
elements in the prime cycle, so the weight of such a prime cycle on a /-set is 
2 l (q + q 2 + - ■ ■ + q l ~ 1 ). Hence the generating function of d^(q) follows from the 
generating function of P-decompositions of type (l fc , h, . . . , l m ), k > 0, U > 2, 
as given by 



\ H) • • • j 'm /„■_-, 



h+- + lr, 



Summing over li, ... ,l m > 2 and k > 0, m > 0, we obtain the right hand 
side of IKTli . I 
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4 A Recurrence Relation 



In this section, we will use the recurrence relation for Eulerian polynomials 
of type B to derive a recurrence relation for the g-derangement polynomials 
dn(q). Applying a theorem of Zhang [15], we deduce that d^(q) form a 
Sturm sequence, that is, d^(q) has only real roots and separated by the 
roots of Moreover, from the initial values, one sees that d^(q) has 

only non-positive real roots. Consequently, d^(q) is log-concave. Although 
the polynomials d„(q) are not symmetric, we show that they are almost 
symmetric in the sense that the coefficients have the spiral property. 

The following recurrence formula (14.11) for B n (q) is a special case of The- 
orem 3.4 in Brenti [3], see also Chow and Gessel [6], which will play a key 
role in establishing a recurrence relation for d^(q). 

Theorem 4.1 We have 

B n {q) = ((2n - l)q + l)B„_i(?) + 2q{l - q)B' n _ 1 (q), n > 1, (4.1) 
where Bo(q) = 1. 

Theorem 4.2 For n > 2, we have 

d*{q) = (2n - l)gd*_i(?) + 2g(l - 9)^(9) + 2(n - l)gC 2 (?)- ( 4 - 2 ) 
Proo/. By (J3J| and (Oil , we obtain 

<(?) = £(-i)"- fe (") 



fc=0 



+ E(-l) n " fc f" _ ]) (((2fc - l)g + l)B fe -i(<z) + 2g(l - 



fc=i 



-gd^fe) + 2q^(-l)«- & ( - f n X ) ) fc5 fc _!( g ) + 25(1 - g^fe) 



fc=i 



-^_ 1 (5) + 2ng^(-ir- fc ^_^5 fc _ 1 (5) 



fc=i 



+ 2g(n - 1) (" _ ^ B k -i(q) + 2g(l - g^fe) 

fc=i ^ ' ' 

= (2n - l)gCi(?) + 2(n - l)<?<£_ 2 (g) + 2g(l - 

as desired. I 

Equating coefficients on both sides of (14. 2p . we are led to the following 
recurrence relation for the numbers d n ^. 
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Corollary 4.3 For n > 2 and k > 1, we /iai>e 

rf njfc = 2kd n _ l ^ k + (2n - 2A; + l)d n _ 1>fc _! + 2(n - l)d n _ 2) fc_i. (4.3) 

From the above relation (14.31) . it follows that d n \ = 2 n for n > 1. The 
recurrence relation (14.21) on <i^(g) enables us to show that the polynomials 
d^(q) form a Sturm sequence. The proof turns out to be an application of 
the following theorem of Zhang [T5] . 

Theorem 4.4 Let f n (q) be a polynomial of degree n with nonnegative real 
coefficients satisfying the following conditions: 

(1) Forn > 2, f n (q) = a n qf n ^ l (q) + b n q(l + c n q)f' n _ l {q) + d n qf n _ 2 (q), where 
a n >0,b n > 0, c n <0,d n > 0; 

(2) For n > 1, zero is a simple root of f n (q); 

(3) fo(q) = e, fi(q) = e^q and fi{.q) has two real roots, where e > and 
ei > 0. 

Then the polynomial f n (q) has n distinct real roots, separated by the roots of 
fn-i(q), n>2. 

It can be easily verified that the recurrence relation (14. 2 p satisfies the 
conditions in the above theorem. Thus we reach the following assertion. 

Theorem 4.5 The polynomials d^(q) form a Sturm sequence, that is, d^(q) 
has n distinct non-positive real roots, separated by the roots of d^_i(q). 

As a direct consequence of the above theorem, we see that the coefficients 
of d^(q) are log-concave. Although the coefficients are not symmetric as in 
the type A case, we will show that they are almost symmetric in the sense 
that they satisfy the spiral property. The spiral property was first observed 
by Zhang [16] in his proof of a conjecture of Chen and Rota [I]. 

Theorem 4.6 The polynomials d^(q) have the spiral property. Precisely, 
forn >2, if n is even, then 

dn, n < d n> i 

and if n is odd, then 

d n , n < d n , 1 < d n , n-l < d n ,2 < d n , n-2 < ' ' ' < ^n, 2+1 < « n , 2=1 < "n, 2+1 ■ 
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Proof. Let 

f(n) 



| — 1, if n is even, 
if n is odd. 



In this notation, the spiral property can be described by the following in- 
equalities 

<d ntk < d n ^ k (4.4) 
for any 1 < k < f{n), and in addition, the inequality 

ci„,f+i<4,f (4.5) 

when n is even. 

We proceed to prove the relations (14. 4p and (14. 5 p by induction on n. It 
is easily seen that (I4.4p and (14. 5 p hold for n = 2 and n = 3. We now assume 
that they hold for all integers up to n. We now aim to show that 

(4.6) 

for any 1 < k < f(n + 1), and it is also necessary to show that when n + 1 
is even, 

d n+ i,2i±3 < d n+h n±i. (4.7) 

For k = 1, we have d n+hn+x -d n+ i A = l-2 n+1 < 0. For 2 < k < /(n+1), 
by the recurrence relation (14. 3p for d n _ k , we have 

dn+l.n+2-k = 2(n + 2 — k)d n , n+ 2- k + (2k — l)rf n , n+l-A: + 2nd n -\, n +l-k, 

(4.8) 

d n+ i >fc = 2kd n<k + (2n - 2k + 3)d nife _i + 2nd n -i >k -i, (4.9) 
d n+ i in+ i_ fc = 2(n + 1 - k)d n , n+1 _ k + (2/c + l)cf n ,n-fc + 2nd n _i in „ fc . (4.10) 
It follows from (j48]) and flUD that 

^n+l,n+2-fc — d n+ i )k = (2n — 2/c + 3)((i n)n+ 2-fc — ^n,fc-l) + 2k(d n ^ n+ \- k — d Utk ) 

+ 2n((i n _i in+ i_fc — C? n _i 5 fc_i) + (d n)n+ 2-k — d ntn+ i- k ). 

By the inductive hypothesis, we see that the difference in every parenthesis 
in the above expression is negative. This implies that for 2 < k < f(n + 1) 

d n +l,n+2-h — d n+ i tk < 0. (4-11) 
Similarly, for 2 < k < f(n + 1), in view of (14. 9 p and (14.101) we find 

d n +i,k — ^n+i,n+i-fc — (2A; + l)(d njk — d n n _ k ) + 2n(d n _i tk _i — d n ^i tn _ k ) 
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+ [2n + 3 — 2k)(d nt k~l — <2 n ,n+l-A;) + (d n ,n+l-k ~ dn,k)- 
Again, by the inductive hypothesis, we deduce that for 2 < k < f(n + 1), 

d n +i,k — d n+ i tn+ i-k < 0. (4-12) 

Combining fl4TTT|) and (TCT2]) gives (|4J| for 1 < k < f(n + 1). 

It remains to verify (14.71) when n + 1 is even. By the recurrence relation 
(14.31) of d nt k, we have 

d n+ i t n±3 = (n + 3)d n ^ n±3 + nd n ^ «±i + 2nd n _ x ^ «±i , 

o? n+ i, »±i = (n + l)d ni n±i + (n + 2)d nt n_i + 2nd n _ h ^ . 

This yields 

rf n+ i,2±3 - d n+h n±i =(n + 2)(d n n ±a - d n ^) + (d n> n±3 - d n>2 ±i) 

+ 2nK_i*-Ci,ti)' 

Again, by the inductive hypothesis, we immediately obtain (14.71) . This com- 
pletes the proof. I 

5 The Limiting Distribution 

In this section, we show that the limiting distribution of the coefficients of 
dn(q) is normal. The type A case has been studied by Clark [7j. It has been 
shown that the limiting distribution of the coefficients of d n (q) is normal. Let 
£ n be the number of type B excedances in a random type B derangement on 
[n]. We first compute the expectation and the variance of £ n . Then we use 
Lyapunov's theorem to deduce that £ n is asymptotically normal. 

Theorem 5.1 We have 

Varf, 

Proof. By the recurrence relation (14.11) for B n (x), we have for n > 1, 
B'Jx) = {2n-l)B n _ l (x) + (2nx-5x + 3)B' n _ 1 (x) + 2x(l-x)B , ;_ 1 (x). (5.3) 



n I ... 
2 + 4 + 0(1) ' 

--- + o(l) 
12 16 v ' 



(5.1) 
(5.2) 
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Since B n (l) = 2 n n\ for n > 0, setting x = 1 in (15.31) gives the following 
recurrence relation for B' n (l): 

B' n (l) = (2n - l)(n - 1)!2"" 1 + (2n - 2)5;_ 1 (1). 

It can be verified that for n > 1, 

*<U = (5-4) 

Moreover, by (15. 3p we get 

B'^x) = (An-Q)B , n _ 1 (x) + (2nx-9x + 5)B , r ' l _ 1 (x) + 2x(l-x)B , r ' l ^ 1 (x). (5.5) 
Setting x = 1 in (15.51) and using (15.4)) . we obtain 

= (2n - 3)(n - l)2"- l (n - 1)! + (2n - 4)^(1). (5.6) 

Then for n > 2, 

B „ (1) = (5 . 7) 

Since 5 n (l) = 2 n n\, in view of the formula (13.31) . we see that 

^(i)=£(-ir fe ("W)- (5.8) 

fc=0 ^ ' 

Let 

n 
fe=0 

Then d^(l) can be written as 2 n n!s n . 

Applying the formula (13 .3p again and using the evaluation (I5.4[) for B' n (l), 
we find that for n > 1 



<'(i) = E(-ir m (^) ■ 2m ~ 1 

m=l ^ ' 



m • m! 



n-1 

^ m!2 m+1 

m=0 



2 n n\ 



n—1 1 n—1 -, 



E(-i) m ^ + iE(-D' 



2 ^ v ' m!2 m 4^ v y 2 m " 1 (m-l)! 

m=0 m=0 



2 n n! / 1 

nsn-i + -s„_ 2 ) • (5.9) 
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Differentiating (13. 3p twice and invoking (15. 7p . we deduce that for n > 2 

, 3m 2 — 5m + 1 



rff(i)=£(-ir m (^V 

m=2 ^ ' 



12 



2"u! ^ ( _ 1)t ,..,.,3m 2 - •')/// I 



12 ^ v ' (n-m)!2 n - 

m=2 v ' 

n-2 „, , 2 



2 n n\ 3(n - m) - 5(n - m) + 1 

19 ' ml?™ 



12 ^ V 7 m!2 r 

Z^l ( y f ir 3n a -5" + l m ~6n + 5 

12 l^ 1 J m\2 m ^ ' (m-l)!2 r 

\m=0 m=0 v y 

n-2 o 

m 



+3E(-i: 



m!2 r * 

m=0 

2 n n! / 1 3 3 

( (3n 2 - 5n + l)s n _ 2 + -(6ra - 5)s n _ 3 + -s n _ 4 - -s n - 

2 n n\ 
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(3n - 5n+ l)s n _ 2 + (3n - 4)s„_ 3 + -s„_ 4 . (5.10) 



It is easy to see that s n _ r /s n = 1 + o(l) for r = 1,2,3,4. From (15. 8p . 
( 15.91) and (15.101) . we conclude that 

Va^ = + E£ n - (E£ n ) 2 = - - - + o(l), (5.12) 

as desired. I 

Given the formulas for the expectation and variance of £„, we will use 
Lyapunov's theorem [T(H Section 1.2] to derive that the limiting distribu- 
tion of £ n is normal. Recall that a triangular array of independent random 
variables k — 1,2, ... ,n, n = 1,2,..., is called a Poisson sequence if 

P{£nk = 1} = Pk, P{ink = 0} = q k , 

where p k — Pk{ n ), Qk — Qk{ n ) and Pk + Qk — 1- Then Lyapunov's theorem 
can be used to derive asymptotically normal distributions [TUl Section 1.2]. 

Theorem 5.2 (Lyapunov) Let 

n n 

V n = ^PkQk, Vn = V' 1 J2(£nk ~ Pk)- 

k=l fc=l 
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If Vn — > oo as n — > oo, t/ien the sequence {f] n } is asymptotically standard 
normal. 

The above theorem enables us to derive the asymptotic distribution of 
the random variable r\ n . 

Theorem 5.3 The distribution of the random variable 

VVar^ n 

converges to the standard normal distribution as n — ► oo. 

Proof. Since the polynomials d^(q) have distinct, real and non-positive roots, 
we may express d^(q) as 

dn (?) = <?(<? + «i)0? + a 2 )---{q + a n -i)- 

Hence the random variable £ n , namely, the number of type 5 excedances 
in a random type B derangement on [n], can be represented as the sum of 
independent random variables 

where £„i, £ n2 , . . . , £ n ,n-i, £n,n form a Poisson sequence with p^ — i+a k ; ^ r 
fc = 1,2, . . . ,n — 1 and p n = 1. On the other hand, the formula ( 15.121) implies 
that Var£ n — > oo as n — > oo. Therefore, from Theorem 15.21 we deduce that 
r] n is asymptotically standard normal. I 
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